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ABSTRACT 

We present an analysis comparing the bulk-flow measurements for six recent 
peculiar velocity surveys, namely, ENEAR, SFI, RFGC, SBF and the Mark III 
singles and group catalogs. We study whether the direction of the bulk-flow 
estimates are consistent with each other and construct the full three dimen- 
sional bulk-flow vectors for each survey. We show that although the surveys 
differ in their geometry, galaxy morphologies, distance measures and measure- 
ment errors, their bulk flow vectors are expected to be highly correlated and 
in fact show impressive agreement in all cases. We found a combined weighted 
mean bulk motion of 330 km s _1 ±101 km s" 1 toward I = 234° ± 11° and 
b = 12° ± 9° in a sphere with an effective depth of ~ 4000 km s -1 . 

Key words: Cosmology: Large-scale structure of the Universe; Galaxies: 
Distances and redshifts 



1 INTRODUCTION 

The analysis of peculiar velocity fields of galaxies and clusters is one of the most effective 
ways of probing mass fluctuations on ~ 100 h -1 Mpc scales (h being the Hubble constant in 
units of 100 km s -1 Mpc -1 ). Studies of peculiar velocities can be used to constrain the am- 
plitude of mass power spectrum on scales others than those probed by redshift surveys and 
those sampled by anisotropies in the CMB (e.g. Zaroubi et al.; 2001; Freudling et al. 1999). 
Originally motivated by the invention of distance indicators based on intrinsic relations be- 
tween galaxy observables (Faber & Jackson 1976; Tully Sz Fisher 1977) that are independent 
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2 Sarkar, Feldman & Watkins 

of redshift, the field has, in recent years, experienced great progress toward constructing large 
and homogeneous redshift-distance samples of galaxies and clusters. Although the analyses 
of early redshift-distance surveys of spiral galaxies (Aaronson et al. 1982) and of elliptical 
galaxies (e.g., Lynden-Bell at al. 1988) led to the development of several statistical methods 
for analyzing peculiar velocity data (Dressier etal. 1987; Kaiser, 1988; Feldman & Watkins, 
1994; Strauss & Willick 1995; Watkins & Feldman, 1995), these studies were hampered by 
the fact that they were based on relatively small and shallow data sets. 

The present generation of redshift-distance surveys consist of larger and higher-quality 
data sets of both spiral (da Costa et al. 1996; Giovanelli at al. 1997a, 1997b; Haynes et 
al. 1999a, 1999b; Karachentsev et al. 2000) and early-type galaxies (da Costa et al. 2000a, 
2000b). These new samples pave the path toward a possible resolution of many discrepancies 
found in earlier samples; however, some quantitative disagreements persist. The earlier sta- 
tistical comparisons of the peculiar velocity fields derived from D n -<r and Tully-Fisher (TF) 
distances found significant difference between them (e.g., Gorski et al. 1989; Tormen at al. 
1993). 

Based on the work of Kaiser (1988), Feldman & Watkins (1994) formulated a linear 
analysis to calculate the theoretical expectation for bulk flow in large-scale surveys as a 
function of the geometry of the survey, its clustering properties, and the assumed power 
spectrum; and applied it to a volume-limited complete sample of 119 Abell clusters (Lauer 
& Postman 1994, hereafter LP) to show that the power spectra considered were inconsistent 
with the LP measurement of bulk flow at the 95%-97% confidence level. The formalism 
was later applied (Watkins & Feldman, 1995) to calculate a measure of correlation between 
results obtained from the Riess, Press, & Kirshner (1995a, 1995b, hereafter RPK) and the 
LP samples. They found that the apparent lack of agreement between the two measurements 
could be explained by the fact that both the LP and RPK samples were dominated by noise 
and incomplete cancellation of small scale motions. 

More recently, in an analysis of the ENEAR sample (da Costa at al. 2000a, 2000b), the 
results obtained by Borgani et al. (2000) pointed toward a statistical concordance of the 
velocity fields traced by spiral and elliptical galaxies, with galaxy distances estimated using 
TF and D n -<r distance indicators, respectively. Following the method described in Feldman 
& Watkins (1994) and Watkins & Feldman (1995), Hudson et al. (2000) also showed that 
the bulk flows measured from four different surveys (SMAC, SC, LPlOk and SNIa) were 
consistent with each other. Reconstruction techniques which compared the velocity to those 

© 0000 RAS, MNRAS 000, 000-000 



Consistency of Bulk Flow Surveys 3 

predicted from galaxy redshift distributions ((Zaroubi et al. 2002; Pike & Hudson, 2005; 
Park & Park, 2006)) show that we can get consistent results for (3 = f2°; 6 : /b where Q m and b 
are the mass density and linear bias parameters, respectively. These papers do not directly 
compare the surveys, but rather calculates best fit for (3 given a velocity field. Pairwise 
velocity comparisons between various samples ((Feldman et al, 2003)) also show that different 
velocity samples produce consistent statistical results for Vt rn and a$ the standard deviation 
of density fluctuations on a scale of 8/i _1 Mpc. The fact that they give consistent results does 
not necessarily indicate consistency between the surveys, since the agreement is indirect and 
there is no attempt to quantify the consistency between the surveys and the field. Agreement 
of a specific statistical characteristic between data sets need not mandate consistent data 
sets. 

In the present paper, we calculate the theoretically expected correlation between the 
estimates of the bulk flows of samples of galaxies in four recent surveys, namely, ENEAR 
(da Costa at al. 2000a, 2000b), SFI (Giovanelli et al. 1994; da Costa et al. 1996), RFGC 
(Karachentsev et al. 2000), SBF (Tonry et al. 2001) and the Mark III catalogs (Willick et al. 
1997). We also introduce an analytical method to calculate the likelihood that two surveys 
both sample the same large scale flows, that is, we study whether measurement errors and 
differences in the distribution and morphology of galaxies in the surveys can statistically 
account for the differences in the directions of bulk-flow vectors. Further, we construct the 
3-dimensional bulk flow vectors for all the surveys mentioned above and calculate the actual 
dot products of the estimates of the bulk flows obtained from these surveys in order to discuss 
their consistency. In § 2, we describe the theoretical background of velocity fields; we explain 
in detail the formulation of our analysis in § 3. A description of the surveys considered in 
our analysis is given in § 4; We then discuss our results in § 5 and conclude in § 6. We argue, 
in the context of the EBW (Efstathiou, Bond & White 1992) Standard cold-dark-matter 
(CDM) power spectrum, that our results also show a consistent statistical concordance. 

2 PHYSICS OF VELOCITY FIELDS 

In the context of the gravitational instability model of structure formation, the motions of 
galaxies are directly related to mass-density fluctuations. On scales of the surveys, the mea- 
sured velocity of galaxies deviate from the Hubble expansion due to local mass distribution. 
Thus peculiar velocity surveys provide a unique method to probe the distribution of mass in 
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the local universe. On scales that are small compared to the Hubble radius, galaxy motions 
are manifest in deviations from the idealized isotropic cosmological expansion 

cz = H r + r • [v(r) - v(0)] (1) 

where c is the speed of light, z is the redshift, Hq is the Hubble constant, r is the distance of 
a galaxy at the redshift z, f is the the unit vector toward the galaxy, and v(r) is the proper 
motion of the galaxy (at position r) with respect to the comoving frame. This component of 
the overall motion of the galaxy is known as its peculiar velocity, arising from the gravitational 
attractions of surrounding overdensities. In Eq. (1), v(0) is the peculiar velocity of the 
observer; It is standard practice to omit this term from the equation and to assume that 
redshift has been corrected to account for the motion of the observer. 

The redshift-distance samples, obtained from peculiar velocity surveys, allow us to de- 
termine the radial (i.e., line-of-sight) component of the peculiar velocity of each galaxy: 

v(r) = r ■ v(r) = cz — H r (2) 

We assume that galaxies trace the large-scale linear velocity field v(r) which is described 
by a Gaussian random field that is completely defined, in Fourier space, by its velocity 
power spectrum P v (k). In the statistical model for peculiar velocities we define the Fourier 
Transform of the line-of-sight velocity r • v(r) such that: 

r-v(r) = — ^ [d 3 kr -kt;(k)e ik - r (3) 

{27r) ,i J 

Due to the isotropy assumed in the Cosmological Principle, the statistical properties of 
kv (k) are independent of the direction of k, and so we may define the velocity power spectrum 

(v(k)v*(k')) = {2ixfP v {k)5 D {k - k'), (4) 

where 5r> is a Dirac delta function, and the averaging on the left-hand-side is over directions 
of k. 

In linear theory, the velocity power spectrum is related to the density power spectrum, 
P(k), by 

Pv(k) = f 2 (n mfl ,n A ) p(k) . (5) 

f(fl m fi, ^a) is the rate of growth of the perturbations at the present epoch and can be 
approximated as (e.g., Lahav et al. 1991): 

/(n mi0 ,n A )«nSJ ) (6) 
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where £l m ,o is the cosmological density parameter for matter at the present epoch. 

The power spectrum provides a complete statistical description of the linear peculiar 
velocity field. It should be noted that the above expressions are valid only on scales sufficiently 
large so that non-linearity can be neglected. In the present analysis, we consider the EBW 
parameterization of the linear CDM power spectrum (Efstathiou, Bond & White 1992) 



where T parameterizes the "shape" of the power spectrum and the overall normalization 
is determined by a s , the standard deviation of density fluctuations on a scale of 8/i _1 Mpc. 
The constant C is determined by the direct relation between <r 8 and the power spectrum. 
For models where the total density parameter Q = 1, the shape parameter is related to the 
density of matter, T = Q m ,oh. In the present analysis we use os = 0.9, T = 0.21 and h = 0.7. 

3 MODELING THE OBSERVATIONAL DATA 

A catalog of peculiar velocities consists of a set of galaxies, labeled by an index n, for 
which we are given positions r n and estimates of the line-of-sight peculiar velocities S n with 
uncertainties a n . For simplicity, we will make the assumption that observational errors are 
Gaussian distributed. Since linear theory only applies on scales comparable to the survey 
size, we focus our attention on the lowest order moments of a Taylor expansion of the 
velocity field v(r). Following Kaiser (1988), we model the velocity field as a uniform streaming 
motion, or bulk flow, denoted by U, about which are random motions drawn from a Gaussian 
distribution with a 1-D velocity dispersion a*. Although this model ignores the fact that 
small-scale motions, including those that are nonlinear, are correlated, it is reasonable to 
assume that they effectively average out on the scales we are considering. Since the value of 
a* is not well determined by linear theory, we will treat it as a parameter with a fixed value 
of 300 km s -1 . We have checked that our results are fairly insensitive to the exact value 
chosen for this parameter. Given these assumptions, the likelihood function for the bulk flow 
components is 



where here and in subsequent equations repeated indices are summed over. The maximum 
likelihood solution for the ith component of the bulk flow is given by 



P(k) = a 2 8 Ck(l + [6.4(fc/r) + 3(A;/r) 1 - 5 + (1.7k/T) 2 ] 1A3 ) 



2/1.13 



(7) 




(8) 



n ° 'i 




(9) 
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6 Sarkar, Feldman & Watkins 
where 



^> = E^ do) 

n u n ' * 

Thus Uj is the cross-correlation between the estimated line-of-sight velocity of the n-th galaxy 
and its position vector. For the catalogs considered, is nearly diagonal, the off-diagonal 
terms being of order 10% of the diagonal ones. 

In the model we are considering, the measured peculiar velocity of galaxy n is related to 
the velocity field at the position of galaxy n by 

S„ = f n ,iVi(T n ) + e n (11) 

where e n is drawn from a Gaussian with zero mean and variance a\ + of. The fact that e n 
is statistically independent of the velocity allows the theoretical covariance matrix for the 
bulk flow components to be written as (Kaiser, 1988) 

= (UP,) = B$ + B§\ (12) 
where the "noise" term can be shown to be (Kaiser, 1988) 

4 } = A? (13) 

and the "theoretical" term can be written as the convolution of an angle-averaged tensor 
window function with the power spectrum 

Rff = 4vr / dk k 2 P v (k) Wfj(k) (14) 



where 



Our main goal in this paper is to figure out whether the surveys we consider are consistent 
with one another. However, even if two surveys are measuring the same underlying velocity 
field, they will not necessarily give the same bulk flow. This is both due to measurement 
errors in the peculiar velocities and the fact that each survey probes the velocity field in a 
different way. This is most clearly seen by observing that each survey has different window 
functions (see below). In order to get an idea of how much correlation is expected between 
the estimates of the components of the bulk flows XJ A and XJ B of any pair of surveys (A,B) 
for a given power spectrum, we can calculate the correlation matrix (\J a \J bS J for the two 
surveys. This is calculated in a similar manner to the covariance matrix, except that the two 
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Consistency of Bulk Flow Surveys 7 
sums in the window function are now over two different surveys 



X /^T ^-k)exp(ik.(r^-0) (16) 

The correlation matrix can then be used to calculate the normalized expectation value 
for the dot product of XJ A and XJ B (Watkins & Feldman, 1995): 

C = N 1 1 1 _ = (cos^), (17) 

WW) (U B U B )f 2 

where 9 is the angle between U A and U B . C should be close to unity for highly correlated 
vectors, zero for vectors that are completely uncorrelated, and -1 if there is a high degree 
of anti-correlation. It is important to realize that C carries information only about the cor- 
relation of the directions of the bulk flow vectors of the two surveys; however, it provides a 
convenient measure of how well the large scale velocity information contained in two surveys 
agree. 

Given a value of C for two surveys (A,B) calculated using a given power spectrum, we 
can estimate the probability that the bulk flow vectors U A and XJ B will be separated by 
an angle greater than some 9 C . Our strategy is to think of the direction of XJ A as scattering 
about the direction of XJ B in a two-dimensional space where 9 is the radial distance. Thus 
we can take 9 to have a x 2 distribution with two degrees of freedom 

P(9)d9 = !L e - e2,2a2 d9. (18) 
a 2 

The probability of measuring a value for 9 greater than 9 C is then 

P{9 >9 C )= / P{9) dO. (19) 

We can estimate the value of a by using the fact that C = (cos 9) ml — \ (9 2 ) for small 9. 
Since our P(9) distribution has the property that (9 2 ) = 2a 2 , we can estimate 

a = VT^C . (20) 



This analysis ignores the small anisotropy in the covariance matrices for XJ A and XJ B , but 
should be sufficient for our purposes. 

© 0000 RAS, MNRAS 000, 000-000 



8 Sarkar, Feldman & Watkins 
4 THE SURVEYS 

The formalism described above can be employed to test the consistency of all velocity field 
surveys. In this study, we have considered the following proper distance catalogs: 

1) Spiral Field I-Band (SFI): This is an all-sky survey (Giovanelli et al. 1994; da Costa et 
al. 1996; Giovanelli et al. 1998; Haynes at al. 1999a, 1999b), containing 1104 late-type spiral 
galaxies with I-Band TF distance estimates. It is a angular diameter limited survey and 
covers a volume out to ~ 70 h" 1 Mpc. 

2) Nearby Early-type Galaxy Survey (ENEAR): This is an all-sky survey probing a volume 
out to ~ 70 h _1 Mpc. Although the survey contains data from different sources, it has been 
conducted by a single group and the data was analyzed by a single procedure and reached the 
same completeness level across the sky (da Costa at al. 2000a, 2000b). The sample contains 
702 independent objects early-type elliptical galaxies and groups of galaxies brighter than 
nig = 14.5 with D n -<x measured distances probing volume similar to the SFI survey. 

3) Revised Flat Galaxy Catalog (RFGC): This catalog (Karachentsev et al. 2000) provides 
a list of radial velocities, HI line widths, TF distances, and peculiar velocities of 1327 spiral 
galaxies. This was compiled from observations of flat galaxies from FGC (Karachentsev, 
Karachentseva, & Pernovsky 1993) performed with the 305 m telescope at Arecibo (Giovanelli 
et al. 1997) confined to the zone 0° < 5 < +38° accesible to the telescope. 

4) Surface Brightness Fluctuation (SBF) : This catalog (Tonry et al. 2001) employs I-band 
surface brightness fluctuation method and consists of 269 galaxies (both spiral and elliptical) 
reaching out to ~ 4000 km s _1 , having a characteristic depth ~ 12 h _1 Mpc. 

5) Mark III catalog of singles: 
and 

6) Mark III catalog of groups: These catalogs (Willick et al. 1997) are a compilation of various 
disparate surveys that were recalibrated and compiled to provide some of the first reasonably 
dense and deep peculiar velocity surveys. The Mark III Catalogs provide the observables for 
each object (i.e. redshift, magnitude, velocity width) and inferred distances derived from 
both the forward and inverse TF or D n — a relations. Distances for both individual objects 
and groups are provided. The singles catalog has 2538 galaxies, while the group catalog has 
1124 groups. The total survey depth is over 100 h" 1 Mpc with homogeneous sky coverage to 
~ 30 h" 1 Mpc. 
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Consistency of Bulk Flow Surveys 9 



The most important point to note here is that the above surveys are by and large inde- 
pendent of each other. They use different distance indicators, selection functions and survey 
geometries, and target different morphology of galaxies. Further, since our formalism weighs 
the galaxy contribution by the distance error, we do not expect distant galaxies to contribute 
much to our results and thus, the homogeneous Malmquist bias correction will not change 
our outcome, (see, e.g. (Ferreira et al. f999; Park & Park, 2006)). 

As for the inhomogeneous Malmquist bias (1MB), since redshift distances are not affected, 
it will result in scattering the distances away from overdensities, leading to an appearance of 
an infall. This effect should not affect the bulk flow, though it may contribute to higher order 
moments. If we had found disagreement between the surveys, it would be very important to 
estimate the magnitude of the 1MB since it may be the cause of the disagreement. Since we 
found agreement, then it is likely that the 1MB is smaller than our errors, since it was not 
sufficient to cause disagreement even though one would expect each survey to have a different 
bias. That said, we have conducted the following experiment: We have cut all galaxies above 
a certain radius (for radii r = 40, 50, 60, 70 Yc 1 Mpc) and performed the analysis described 
above. Removing far away galaxies did not change the results significantly and the direction 
of the flow for these radii was within a standard deviation of the results quoted in this paper. 
These results are expected, since the galaxies are distance-error-weighted and thus, galaxies 
close to the edge of the surveys contribute relatively little to the overall flow. We thus confirm 
our expectation that prominent features close to the edge of the surveys (e.g. Pieces-Perseus 
region) do not produce a spurious 1MB bulk motion and systematically bias our results. 



5 RESULTS 

We now present the estimates of the actual bulk flow vectors for all the different surveys. 
We use Eq. (9) to construct the Cartesian components of the full three-dimensional bulk 
flow vectors. Our results are tabulated in Table 1. The uncertainties given for the bulk flow 
values in Table I are obtained from the noise part of the covariance matrix, (Eq. 13). 
Since is nearly diagonal for all of the surveys, we took the uncertainties in the Ui to be 



approximately independent, so that the uncertainties are taken to be ou i = \J R u . These 
uncertainties are dominated by measurement errors in the individual velocities, although 
they also have small contributions from cr* and cosmic scatter. 
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Table 1 



Survpv 


Method 


N 


Effective Depth 




U„ 

V 










(km s _1 ) 


(km s _1 ) 


(km s _1 ) 


(km s _1 ) 


SFI 


TF 


1104 


~ 4000 


43± 37 


-145± 35 


57± 26 


ENEAR 


D n -a 


702 


~ 4000 


154± 50 


-246± 44 


17± 39 


RFGC 


TF 


1280 


~ 6000 


235± 38 


-232± 39 


31± 29 


SBF 


SBF 


280 


~ 2000 


249± 58 


-262± 46 


163± 29 


Mark III singles 


Various 


2538 


~ 4500 


202± 26 


-230± 24 


31± 22 


Mark III groups 


Various 


1124 


~ 4500 


247± 34 


-386± 31 


79± 26 



The Cartesian components of the full three-dimensional bulk flow vectors as well as the distance estimator method, 



the number of data points and the effective depth for each survey. 

Each catalog surveys a different volume of space, it samples a small subset of the un- 
derlying population of galaxies, and uses independent spatial sampling techniques. Since the 
universe is homogeneous and isotropic on large scales and the sample volumes of the various 
data sets strongly overlap, these surveys react to the same underlying large-scale mass distri- 
bution. However, the contribution from small scale nonlinearities differs from one catalog to 
another depending on the particulars of the survey. The window functions of these surveys, 
therefore, differ from one another, particularly on small scales. As a consequence, we don't 
expect the measured bulk flows from these surveys to be identical, even in the absence of 
peculiar velocity errors. Thus the bulk flow components listed in Table 1 are not strictly 
comparable. However, if the samples are true representations of the flows of the underlying 
populations, they should be correlated. 

For illustration purposes (Fig. 1) we show the tensor window functions (Eq. 15) for each 
of the surveys. We chose to show the W^ x in the k X) k y plane (we have not performed the 
angle-averaging that is given in the equation). The figure illustrates the fact that although all 
of the surveys have a similar central peak around k = which samples the large-scale power 
in a similar way, each survey samples the region of larger k differently. These differences in 
the window functions tend to decrease the correlation between the bulk flow vectors of the 
surveys. 
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SFI MarkSsmgles 




Figure 1. The logarithm of the trace of the normalized square tensor window functions in the k x — k y plane from the six 
surveys we used in this study. As we can see, the central peak is similar in all surveys, suggesting that the large scale (small 
k) power is probed in a similar fashion. However, as we move away from the center, each survey samples the underlying power 
differently. 



In Table 2, we show the value of C = (cos#) (Eq. 17) for each pair of surveys (A,B), 
together with the inferred a (Eq. 20) for the probability distribution of 9. We can see from 
these values that the directions of the bulk flow vectors for all of the surveys are highly 
correlated. We also show the angle 9 C between the measured U A and XJ B , and the probability 
of measuring an angle this large or larger, P(9 > 9 C ) (Eqs. 18-19). These show that in general 
the results are consistent with one another for all pairs. 
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Table 2 



Survey 


cos(#) 


(cos#) 


a 


9c 


P(9 > 9 C ) 


SFI-ENEAR 


0.92 


0.91 


0.30 


0.40 


0.42 


SFI-RFGC 


0.85 


0.90 


0.31 


0.55 


0.21 


SFI-SBF 


0.91 


0.78 


0.47 


0.44 


0.64 


SFI-Mark III s 


0.88 


0.89 


0.33 


0.49 


0.33 


SFI-Mark III g 


0.95 


0.92 


0.27 


0.33 


0.49 


ENEAR-RFGC 


0.97 


0.86 


0.37 


0.23 


0.82 


ENEAR-SBF 


0.92 


0.80 


0.45 


0.41 


0.66 


ENEAR-Mark III s 


0.99 


0.88 


0.34 


0.17 


0.89 


ENEAR-Mark III g 


0.99 


0.92 


0.29 


0.11 


0.93 


RFGC-SBF 


0.95 


0.77 


0.48 


0.33 


0.79 


RFGC-Mark III s 


1.00 


0.82 


0.42 


0.07 


0.99 


RFGC-Mark III g 


0.97 


0.86 


0.38 


0.23 


0.83 


SBF-Mark III s 


0.95 


0.68 


0.56 


0.32 


0.85 


SBF-Mark III g 


0.95 


0.82 


0.42 


0.31 


0.76 


Mark III s-Mark III g 


0.99 


0.92 


0.28 


0.17 


0.83 



For each pair of surveys we show the value of the cosine of the angle (9) between their bulk flows, 
the expectation value of their dot product C, the inferred width a for the probability distribution 
of 9, the critical angle 9 C and the probability of measuring an angle greater than theta c . 



To test our theoretical results and see in more detail the exact distribution of the bulk flow 
vectors, we conducted numerical experiments with the data. In one experiment we perturbed 
galaxies' positions, and hence also the peculiar velocities, using the reported errors in the 
distance measurements. Essentially we took each catalog and performed 1,000 Monte-Carlo 
realizations of the data using the measurement errors as the width of a Gaussian about the 
mean distance - the proper distance reported. In another experiment we used the diagonal 
elements of the "noise" part of the covariance matrix as the variance of the individual 
components of the bulk flow vectors and did another 1,000 Monte-Carlo realizations where 
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we drew values from a Gaussian distribution, N(/z,E 2 ), with the individual components being 
taken as the mean /i and identifying £ 2 as the variance obtained from R^f [refer to Eq. (13)]. 
Both of our methods for calculating the spread in the components of the streaming solution 
yield similar results and agree with the errors reported in Table 1. 

We thus compared the bulk-flow vectors in two distinct ways, one was to measure the 
vectors directly from the data, monte-carlo the results and compare them, the other was 
to use the power spectrum to estimate the probability that two surveys measure different 
directions for the bulk flow. Since we used the power spectrum to calculate the likelihood 
that two surveys both sample the same large scale flows, one has to discuss the effect of 
cosmic scatter since that is part of the variance. However, surveys that have good all-sky 
coverage, as our surveys do, will have nearly identical contributions to their bulk flows from 
large scales; the dominant part of the variance comes from small scale effects such as galaxy 
noise and distance measurement errors. Thus only a very small part of the differences in 
direction of the bulk flows will come from large scale effects such as cosmic scatter. In order 
to remove the cosmic scatter part of the variance, we need to make a prior assumption, that 
we know the local velocity field from other sources, specifically reconstruction of density 
surveys. We chose not to use this prior and we show, in a model independent way, that 
velocity field surveys are consistent with each other and that they do sample the underlying 
large-scale in an unbiased, robust way. 

The bulk flow vector direction for each survey is given in Fig. 2 where we plot the results 
from perturbing the distances, as described above, in the Aitoff-Hammer projection. It is 
clear from the figure that the bulk flow vectors for all surveys cluster about the same direction 
in the sky. Although the bulk flow components are not strictly comparable, it was hard to 
resist combining the results for all six catalogs to get an estimate of the mean bulk flow of 
a sphere with an effective depth of ~ 4000 km s _1 to be approximately 330 km s _1 ±101 
km s -1 toward / = 234° ±11° and b = 12° ± 9° where / and b are the galactic longitude and 
latitude respectively. The value of the combined bulk flow vector was calculated by weighing 
the results by their errors, finding the weighted mean of the values for each survey. We would 
like to emphasize that this result should be taken with a grain of salt, since the bulk flow is 
volume dependent and these surveys strongly overlap but do not strictly occupy the same 
volume. We would also like to point out that the overall agreement between the bulk flow 
vectors of the different surveys may suggest that the internal sheer for the flows should be 
small, a conclusion we are testing in an upcoming paper (Watkins & Feldman, 2006). 
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23456 23456 

D E (xlOOO km s" 1 ) 

Figure 2. The top panel depicts the spread in the individual components of the bulk flow vector for all the six surveys, depicted 
in an Aitoff— Hammer Galactic projection, the cross indicates the weighted mean bulk flow direction. The bottom left and right 
panels show the galactic longitude and latittude, respectively as a function of the estimated effective depth (Dg) of the survey. 
The solid and dashed lines mark the mean and standard deviation, respectively, of the weighted mean results from the six 
catalogs of I = 234° ± 11° and b = 12° ± 9° 

6 CONCLUSION 

We have presented statistical analyses of the bulk flow measurement for six proper distance 
surveys. We have shown that the estimates of bulk flows obtained from these surveys are 
expected to have a high degree of correlation. Further, we have constructed the actual three 
dimensional estimates of the bulk flow vectors and shown that consistent results are ob- 
tainable from independent distance indicators, once they are applied to uniformly selected 
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samples of galaxies. We find no statistically significant differences between the velocity fields 
mapped by different morphologies, galaxy types or distance indicators. 

We would like to stress that one should not be putting too much emphasis on comparing 
the components of the bulk flow vectors directly, since they are not really comparable. This 
is especially true since the error bars reflect only the statistical and measurement errors and 
do not capture the differences in the bulk flows due to the fact that they probe the power 
spectrum differently. Basically, we don't expect the bulk flow components to agree within 
the error bars shown in Table 1. However, when we look at the direction of the bulk flow 
vectors, they do agree with each other remarkably well. Thus we conclude that all bulk flow 
measurements are consistent with each other given the errors as long as we allow for small 
scale alliasing and incomplete cancelations. A rough estimate of the (weighted mean) bulk 
flow from all surveys gives a flow of 330 km s _1 ± 101 km s" 1 toward I = 234° ± 15° and 
b = 12° ±9°. 

This study clearly supports the notion that we have reached an era where velocity field 
data is consistent and robust across morphological types, selection criteria, survey geometry 
etc. Results from independent catalogs probe the same underlying large-scale power, though 
are subjected to different small-scale fluctuations. Unlike earlier, sparser surveys, the newer 
proper-distance surveys provide us with a dynamical probe of the large-scale structure 
which we can add to our growing arsenal of data with confidence that the results reflect the 
cosmology we probe. 
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